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ABSTRACT 

We calculate the production cross-sections for the vector and pseudoscalar bound 
states of two gluinos. It is shown that existing and future colliders imply a realistic 
chance of observing gluinonium as a narrow peak in the two-jet invariant mass 
spectrum. With an integrated luminosity of 0.2 fb~^ at the Tevatron, and the 
high efficiency for tagging heavy quark jets at CDF, one should be able to detect 
vector gluinonium for gluino masses up to about 170 GeV; or up to about 260 
GeV for an upgraded Tevatron with a centre of mass energy of 2 TeV and an 
integrated luminosity of 1 fb^^. The significantly higher energy and luminosity of 
LHC should allow pseudoscalar gluinonium to be detected for gluino masses up to 
about 1500 GeV for an assumed luminosity of 200 fb^^. These results are insensitive 
to the details of supersymmetry models, provided that R-parity is conserved and 
the gluinos are lighter than the squarks. In addition, gluinonium detection implies 
a relatively accurate measure of the gluino mass, which is difficult to determine by 
other means. 



1 Introduction 



In most supersymmetry models |l|], the ease with which gluinos can be detected 
depends critically on the relative masses of gluinos and squarks. If the gluino mass 
rrig is greater than the squark mass rriq, the gluino can decay to a squark and 
antiquark (or an antisquark and quark). This is a strong decay with a decay width 
of order a^mg, which has prominent signatures, so that the Particle Data Group 
give a current experimental lower limit of rrig > 218 GeV. On the other hand, if the 
squarks are heavier than the gluino and R-parity is conserved, the gluino can only 
decay weakly, as shown in Figure 1, with a width of order [|^ 

^ A8tt \mq 

The gluino is then more difficult to detect and the mass limit is reduced to 100 
GeV 0. However, in the latter case there is the possibility of observing the gluino 
indirectly, by detecting a gluino-gluino bound state. This has the advantage that the 
conclusions which can be drawn from a search for such states hold in a very wide 
class of supersymmetry models. In addition, the detection of such a state would 
lead to a relatively precise determination of the gluino mass, which could not be 
obtained easily by observing the decay products of the gluino itself, as some of these 
escape undetected. Here we investigate the possibility of detecting such states as 
narrow peaks in the di-jet invariant mass distributions, in the light of the significant 
improvements in heavy quark jet tagging and beam luminosity which have been 
achieved at the Tevatron, and in the light of the high energies and luminosities 
expected at the LHC. 

2 Gluinonium 

Gluinos are heavy spin- 1/2 fermions which interact strongly via gluon exchange. 
Like gluons, they have no basic electroweak couplings. Hence, in contrast to top 
quarks, they live long enough to form two gluino bound states provided rrig < rriq 
as indicated above. The resulting spectrum of states is called gluinonium and has 
been studied in a number of papers Q-Q]. For heavy gluinos, the low lying states 
are essentially Coulombic, with masses M ^ 2m^, and they have much in common 
with heavy quarkonia. However, there are two important differences. 
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Firstly, the gluinos form a colour octet, so that there are six possible colour 
states for a two gluino system, corresponding to the well-known decomposition 

8®8 = 1©8s©8a©10©T0©27. (2) 

The single gluon exchange potential V{r) is attractive for the first three, and is 
related to the corresponding potential Vq{r) for qq pairs by 

/9 9 9\ 

for the (1,8s, 8a) states respectively. 

Secondly, the gluino is a Majorana fermion with just two degrees of freedom. It 
is its own antiparticle and the Pauli principle applies to the two gluino system, thus 
forbidding certain states. 

The resulting spectra of low lying states are shown in Table 1. The allowed 

colour singlet 1 and symmetric octet 85 states have the same values as the char- 
monium states with C = 1; while the allowed antisymmetric colour octet 8^ states 
have the same values as the charmonium states with C = —1. In particular, the 
lowest lying colour singlet and symmetric colour octet states are the pseudoscalars 
?7? and 77? with = 0~, while the lowest lying antisymmetric colour octet state is 
vector gluinonium tp^ with = 1~. 

The properties of these low lying states have been discussed previously 0-0. 
Here we simply summarize the main results. All three L = states decay pre- 
dominantly via gluino-gluino annihilation as shown in Figure 2. Specifically, the 
pseudoscalars r/?'^ decay mainly to two gluons ^ Q with decay widths 

while vector gluinonium ip^ decays predominantly into qq pairs [0 with a decay 
width 

- qq) = N.'^lRsm' ^ iV,ga|M . (6) 

Here M ^ 2mg is the gluinonium mass, Ng is the number of quark flavours for which 
the decay is allowed kinematically, and quark masses have been neglected. i?o,8(0) 
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are the radial wavef unctions of the corresponding two gluino systems evaluated at the 
origin. For simplicity, they have been approximated above by their values for a pure 
Coulombic potential, which is a good approximation for very heavy gluino masses. 
However, in our actual calculations we have also incorporated a linear confinement 
term |], which increases the predicted widths by 40% (20%) for M=200 GeV (600 
GeV). Although much larger than the decay width (0) of the constituent gluinos, 
these decay widths are still very small compared to the gluinonium mass M 2mg. 
The size of all three states is of order = A{asM)^^, which is much smaller 
than the confinement length, thus justifying the relative stability of the colour octet 
states(see i)- 

To summarise, vector gluinonium is a heavy compact object which behaves 
rather like a heavy gluon, except that its coupling to quarks is much stronger than 
its coupling to gluons. Hence it is most readily produced via qq annihilation, and 
the Tevatron is at least a potentially promising place to look, being a source of both 
valence quarks and valence antiquarks. In contrast, the pseudoscalar states rj^ and 
r^l couple predominantly to gluons, and can be produced equally well in both pp and 
pp collisions via the gluon-gluon fusion mechanism. Their production cross-section 
increases more rapidly with energy than that for vector gluinonium, and there is 
more chance of detecting them at the LHC We discuss these possibilities in turn. 

3 Vector Gluinonium at the Tevatron 

In this section we consider the possibility of observing vector gluinonium as a peak 
in the di-jet invariant mass distribution at the Tevatron. This possibility has been 
enormously enhanced since the previous study [0], by improvements in the efficiency 
of heavy quark jet tagging at CDF and by the increase of the Tevatron luminosity. 
The vector gluinonium is produced and decays in pp collisions via the subprocess 

q + q^^l ^q + q, Q + Q (7) 

where we use the symbols q = u,d,s and Q = c,b,t to distinguish light and heavy 
quarksQ 

^Obviously, t-quarks contribute only if the gluino is heavy enough, and even then for the range of 
gluino masses accessible at the Tevatron this contribution is strongly suppressed by the available 
phase space. So, in our rough estimates we consider just two "taggable" heavy quark flavours, 
though in the full simulation correct mass-dependent branching fractions have been used. 
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The nature of the background depends on the value of the ratio M/ ^/s. At 
Tevatron energies, the range of interest hes mainly in large values Mj^ > 0.2, 
where the luminosity of colliding qq pairs prevails over the gluon-gluon pairs. In 
this region, the main sources of two-jet background are the subprocesses 

q + q'^ g + g, q + q, Q + Q (8) 

where the first two have the angular dependence 



OC (1 — COS d 



d cos 9* 

That is, they peak sharply at cos^^* = ±1, where 6* denotes the scattering angle in 
the centre of mass frame of the two jets. In contrast, the signal from the subprocess 

(0) has a much weaker dependence 

d^ 1 , 2/1* 

— ~ 1 + cos 9 . 

a cos u* 

Hence, a cut | cos^^*| < z should improve the signal-to-background ratio. However, 
if the cut z is too low, the signal detection efficiency is low. The optimum value of 
z in this case is close to 2/3. 

The usefulness of heavy quark tagging is clearly brought out by considering the 
production ratios for the various final states in both the signal (|^ and background 

(1) . The relative contribution of the three background subprocesses in (H) at small 
I cos 9* I is given by [0, ^ 

gg:qq:QQ = U:Q5: 6, (9) 

while for the signal (|^) one has 

gg:qq:QQ = 0:3:2. (10) 

Hence by tagging the heavy quark jets one reduces the background by a factor of 
85/6 ~ 14, while retaining 40% of the signal. At larger |cos^*| the gain in the 
signal-to-background ratio caused by heavy quark jet tagging will be even larger. 

At smaller gluinonium masses M ^ 2mg < 200 GeV, initial gluons contribute 
much more significantly to the background, even with heavy quark jet tagging, 
through the subprocess 

g + g9£^Q + Q, (11) 
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This makes the signal-to-background ratio hopelessly small for any reasonable jet-jet 
invariant mass resolution. However, this region of gluino masses is already covered 
by other methods. 

So, most of the two-jet QCD background at large invariant masses arises from 
light quark and gluon jets, and the signal-to-background ratio can be significantly 
enhanced by triggering on heavy quark jets 0. To check that this makes the de- 
tection of vector gluinonium a viable possibility at the Tevatron, we have simulated 
both the gluinonium signal and the 2-jet QCD background using the PYTHIA 5.7 



Monte Carlo event generator [11|. The vector gluinonium production and decay was 
simulated by exploiting the fact that tp^ behaves much like a heavy Z' with axial 
current and lepton couplings set to zero and a known mass-dependent vector cur- 
rent coupling to quarks, chosen to comply with the corresponding decay width after 
taking into account appropriate colour and flavour counting. This effective coupling 
included the non-Coulomb corrections mentioned earlier, and an enhancement due 
to the fact that numerous radial excitations of the which could not be separated 
from it for any reasonable mass resolution, will also contribute. These yield an 
overall factor of between 1.8 and 1.6 depending on M, and the resulting effective 
vector coupling ay falls exponentially from ay = 0.225 at M = 2m^ = 225 GeV to 
ay = 0.172 at M = 2mg = 450 GeV. This signal sits on a much larger background, 
which has been simulated on the assumption that it arises entirely from the leading 
order QCD subprocesses for heavy quark pair production (H) and ([TI|). A constant 
K factor K = 2.0 has been used for both signal and background. 

The cross section for vector gluinonium production at the Tevatron is shown in 
Fig. 3 for y/s = 1.8 TeV. Only decays into heavy quark-antiquark pairs were taken 
into account, and a heavy quark jet tagging efficiency of 50% was assumed. The cut 
on the jet angle 6* in the two-jet c.m.s. frame is | cos^^*| < 2/3, and the cut on jet 
rapidity y is \y\ < 2.0. The analysis shows that with an integrated luminosity of 200 
pb^^ one can hope to see the gluinonium signal from gluinos with masses up to 140 
GeV as a 5 standard deviation peak, and the signal from gluinos with masses up to 
170 GeV as a 3 standard deviations peak. The signal-to-background ratio is around 
7 — 10% at the peak for the assumed two-jet invariant mass resolutions of 25 GeV, 
30 GeV and 38 GeV at M = 225 GeV, 320 GeV and 450 GeV respectively |, |l| . 



For an upgraded Tevatron with the energy increased to 2.0 TeV and an integrated 
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luminosity of 1000 pb~^, the reach is significantly higher. The cross-section obtained 
with the same cuts and the same assumptions about the detector resolution and 
heavy quark jet tagging efficiency is shown in Fig. 4. In this case one can hope to 
see the gluinonium signal from gluinos with masses up to 220 GeV as a 5 standard 
deviation peak, and the signal from gluinos with masses up to 260 GeV as a 3 
standard deviation peak. Finally, we note that the statistical significance of the 
peak is essentially inversely proportional to the two-jet invariant mass resolution, so 
the reach can be significantly extended if some way is found to improve the latter. 

4 Pseudoscalar Gluinonium at the LHC 

Pseudoscalar gluinonium can also be detected as a peak in the the di-jet invariant 
mass distribution. However, both rf^ and 77? decay predominantly into gluon jets, 
which are also the main component of the background. Quark jet tagging does not 
help in this case, and one must rely on high statistics and good resolution to separate 
the signal from the background. These can hopefully be achieved at LHC, as we 
shall now discuss. 

Pseudoscalar gluinonium is produced, as we have noted, via the subprocess 

9 + 9^r,T^9 + 9 (12) 

while the 2-jet background in pp collisions is dominated by the gluon scattering 
QCD subprocess 

9 + 9 — '9 + 9 (13) 
with some contribution from the valence light quark elastic scattering subprocesses 

g + q — *9 + Q, Q + Q — 'Q + q (14) 

at higher values of M/ ^/s. Once again, the angular dependence of the background 
has large peaks in the forward and backward directions, da /d cos 6* oc (1 — cos^ 6*)~'^, 
while the signal is isotropic in the two jet centre of mass frame. Thus a cut | cos^^*| < 
z with optimal z close to 0.6 will increase the signal-to-background ratio. 

To check that pseudoscalar gluinonium detection is a viable possibility at the 
LHC, we have again simulated the QCD background, including all the above sub- 



processes, using the PYTHIA 5.7 generator Gluinonium production has not 
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been explicitly implemented in PYTHIA, but all we need is a pseudoscalar object 
of the right mass which is produced by gluon gluon fusion and decays to two gluons. 
In other words, pseudoscalar gluinonium behaves like a pseudoscalar state of heavy 
quarkonia with appropriately adjusted mass, couplings and width. In this way the 
signal can also be simulated in PYTHIA, where we have again used a K-f actor of 2.0 
and an enhancement factor of 1.5 to allow for the fact that numerous radial excita- 
tions of the 7]^ and t]?, which could not be separated from them for any reasonable 
mass resolution, will also contribute. In addition one can not separate the signals 
from 7]^ and 77? either, so a single effective width P has been used to describe the pro- 
duction of all pseudoscalar gluinonium states. This effective width was calculated 
by summing properly weighted widths and (^, with non-Coulomb corrections 
and higher resonance contributions included, giving P = 0.85 GeV, 1.1 GeV and 1.6 
Gev for M =200 GeV, 600 GeV and 2 TeV, respectively. 

The resulting cross section for pseudoscalar gluinonium production at the LHC 
is shown in Fig. 5. The cut on the jet angle 9* in the two-jet centre of mass frame 
is |cos^^*| < 2/3, and the cut on jet rapidity y is \y\ < 2.5. The estimated two- 
jet invariant mass resolution varies from 22 GeV at M = 300 GeV to 32 GeV at 



M = 600 GeV and 50 Gev at M = 2 Tev |T^. The signal-to-background ratio is 
smaller than in the previous section, of order 1-2%, but the much higher luminosity 
of LHC leads none the less to a higher reach in gluino mass. As can be seen in Fig. 
5, with a luminosity of 10^ pb~^ one expects to reach the gluino masses 800 (1000) 
GeV at 5 (3) standard deviations, while for the highest likely integrated luminosity 
of 2 ■ 10^ pb^^ gluino masses up to 1300 (1500) GeV could be covered. 



5 Conclusions 

We have studied the production characteristics of various bound states of two gluinos 
(gluinonium) for existing and future hadronic colliders. We conclude that gluino- 
nium states can be detected as narrow peaks in the di-jet invariant mass spectra, 
effectively complementing more traditional gluino searches, in the case when the 
gluino is lighter than the squarks. 

In pp collisions one expects copious production of vector gluinonium, which de- 
cays predominantly to qq pairs. The high efficiency of the heavy quark jet tagging 
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at CDF, together with the boost of the Tevatron luminosity, should allow one to 
reach gluino masses of 140-170 GeV at = 1-8 TeV, and 200 pb"^ and 220-260 
GeV at = 2.0 TeV and 1000 pb~^, with realistic efficiencies, resolutions and 
experimental cuts taken into account. 

In pp collisions pseudoscalar gluinonium production is quite large, again allowing 
one to look for a narrow peak in di-jet invariant mass spectra in spite of the high 
background from QCD scattering subprocesses. The high luminosity of the LHC, 
together with the high design resolutions of the ATLAS and CMS detectors, should 
allow one to reach gluino masses of 800 — 1000 GeV at ^/s = 14 TeV and 10^ pb~^, 
while for the highest likely integrated luminosity 2 ■ 10^ pb^^ the reach increases to 
1300 - 1500 GeV. 

Finally, we note that if gluinonium is observed, its mass can be determined with 
an accuracy of a few GeV. This makes gluinonium detection by far the most precise 
way of measuring the gluino mass, since the binding energy inside gluinonium is 
very small and calculable. 
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Table 1. Spin-parities J for the low lying states of gluinonium with L < 2. 
The three columns correspond to the colour singlet state 1 and the symmetric and 
antisymmetric colour octet states 85 and 8^ respectively. 
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Figure 1: The dominant decay mode for gluinos, if they are hghter than the squarks, 
where x° is a neutrahno. 




Figure 2: The dominant decay mechanisms for (a) pseudoscalar and (b) vector 
gluionium. The dominant production processes in hadron coUisions are obtained by 
reading the diagrams from the right. 
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Figure 3: The calculated production cross section of vector gluinonium in pp colli- 
sions at 1.8 TeV. The sohd and broken horizontal hnes indicate the cross sections 
corresponding to a statistical significance at the peak of 5 and 3 standard deviations 
respectively, for a luminosity of 200 pb"^. (See the text for the cuts and resolutions 
used) . 
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Figure 4: The calculated production cross section of vector gluinonium in pp 
collisions at 2.0 TeV. The solid and broken horizontal lines indicate the cross sections 
corresponding to a statistical significance at the peak of 5 and 3 standard deviations 
respectively, for a luminosity of 1 fb~^. (See the text for the cuts and resolutions 
used) . 
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Figure 5: The calculated production cross section of pseudoscalar gluinonium in 
pp interactions at 14 TeV. The sohd and broken horizontal lines indicate the cross 
sections corresponding to a statistical significance at the peak of 5 and 3 standard 
deviations respectively, for two values of the luminosity: 10^ pb~^ and 2 • 10^ pb~^. 
(See the text for the cuts and resolutions used). 
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